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Введение. При построении систем динамической диагностики узлов трения в роторных системах 
типичной задачей является выяснение взаимосвязи изменяющейся геометрии контактирующих 
поверхностей с наблюдаемыми координатами. Для этого, прежде всего, необходимо создать ма- 
тематическую модель динамики диагностируемого объекта. Рассмотрим методику построения мо- 
дели роторной системы. 

В большинстве случаев роторные системы представляют собой тела (валы), имеющие две 
различные точки закрепления в подшипниках. Для математического описания роторной системы, 
совершающей пространственное движение, необходимо задать две системы дифференциальных 
уравнении: систему, определяющую динамику движения центра масс, а также динамику сфериче- 
ского движения. 

Уравнение Эйлера в параметрах Родрига-—Гамильтона. Обозначим неподвижную прямо- 


угольную систему координат как Х | а | ‚ где х, — орты соответствующих осей (рис. 1). 
Кроме того, свяжем с валом подвижную прямоугольную систему координат У = Я У, У где 


у, — орты соответствующих осей. Причем, выберем У так, чтобы начало системы координат сов- 


падало с центром масс вала ©, а оси совпадали с главными осями инерции вала в точке С. То- 
гда, в соответствии с теоремой о движении центра масс [1], первая система уравнений примет 
ВИД: 


тс” =Е^, $=13, (1) 


1 1 


где с“ - координаты точки © (центра масс) в Х; Е” - проекции главного вектора внешних 


сил РЁ на оси в Х !; т - масса вала. 





Рис. 1. Обобщенная роторная система 


44 


Вестник ДГТУ. 2012. № 1(62), вын. 1 


Вторая система дифференциальных уравнений, определяющая сферическое движение 
твердого тела, как известно, описывается уравнениями Эйлера [1]: 


Л + 65. 6 Х =.) = М, 
1,65 + в в; (Л -/,)=М., (2) 
Льфз +565 (Л, -Л)=М,, 


— 


где М, ‚в, - проекции соответственно главного вектора внешних моментов М и угловой час- 


тоты на оси У; Л, — моменты инерции. 


Очевидным недостатком уравнений (2) является то, что они представлены в подвижной 
системе координат. Кроме того, из них явным образом не следуют динамика положения вала, так 


как <, задает лишь его «скорость изменения». Для устранения последнего недостатка необходи- 
мо задать способ представления вала в пространстве, т. е. такую систему переменных, которая 
бы однозначным образом определяла ориентацию вала, и выразить «, через данную совокуп- 


ность переменных. Существуют различные способы решения данной проблемы. В качестве таких 
переменных нами использованы так называемые параметры Родрига-Гамильтона. 

В соответствии с теоремой Эйлера-Даламбера, абсолютно твердое тело из произвольного 
начального положения может быть переведено в любое другое положение (при сферическом 
движении) посредством лишь одного конечного поворота вокруг некоторой оси ОЁ на угол ф. 


Это свойство положено в основу кватернионной кинематики сферического движения твердого 
тела [2]. Исходя из нее любое положение твердого тела (относительно некоторого начального) 


3 
может быть описано посредством кватерниона А=+ > 1А,, где ^, - параметры Родрига- 
1=1 
Гамильтона, 1, — кватернионные единицы. Кватернионные единицы могут быть ассоциированы с 
ортами х,, х,, х, ‚ т. е. будет иметь место х, =1,. При этом говорят, что кватернион представлен 


в базисе Х. В дальнейшем рассмотрены лишь кватернионы в Х. В кватернионной кинематике на 
параметры Родрига-Гамильтона налагают условие нормировки: 


Мо +5 + А: =1. (3) 
Для любого вектора Г, представленного координатами Г, в неподвижном базисе Х, 
можно найти его координатное представление в базисе У по формуле: 
Г, =РИх, (4) 
где Р - матрица направляющих косинусов, определенная как 
№ + — > Л. 2(7.^. + №ойз ) 2(^.^. ЕЛЬ) 
= 2(^.^. Аз ) м -№ +А> Л) 2(^.,^ +А. А ) 
2(7%.^. +. ) 2(%.^. — А ) № -М-^, +^ 
Установим связь между проекциями вектора угловой частоты « на подвижные оси У и 


параметрами Родрига-Г амильтона. Для этого используем кинематическое уравнение Пуассона: 


А-5^зб, (5) 


где символ «о » обозначает операцию кватернионного умножения. 
Используя матричное представление кватернионов [2], можно получить координатную за- 
пись данного выражения: 


в =2 (2%, М = Л.А, +2. А, , 
6% =2 (А. А, +24 Аь — 2,9 ), (6) 
в =2 (№40, -№4А, +А,А, - АА ). 
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Подставив уравнения (6) в уравнения динамики сферического движения (2), получим сис- 
тему уравнений динамики сферического движения в параметрах Родрига-Гамильтона. Причем 
условие нормировки (3) необходимо рассматривать в качестве одного из первых интегралов дви- 
жения. Однако полученная система имеет неудобную форму для проведения численного модели- 
рования. Получим более компактную форму. Из формулы (6), а также из условия нормировки, 


несложно видеть, что имеет место: 
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где 


Со <> 5 


® №5 А+ № 
№ № —^> ” , (7) 
А. м. №, А 
А. = Ло Л: 
№. М сб № 
—. о № —^. А А 
а 
А М № А, в 
^ №5 А Ао 
А АА, 
—^ о № А. | 
оО м № |, 


под символом «® » подразумевается операция поэлементного умножения. 


Операция поэлементного умножения определяется следующим образом: если 2 =| ау. и 


В = ГВ, | ‚ то матрица, определенная как С = А» В , будет состоять из элементов с, =@ 


Ь 


Продифференцировав (7), с учетом (3) получим: 


0 ^ 
& |_,|-^ 
(У | _ 

(> Е 
65 Аз 
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А А ПА. 
 -А Л 
, | —Мь, (9) 
№ № Л, 
—А. Ло А. 


где М, - вспомогательное обозначение, определяемое следующим образом: 


М, =|-2(%+м ++) 000. 


Введем следующие обозначения: 
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Тогда, с учетом вышесказанного, уравнение Эйлера (2) в параметрах Родрига-Гамильтона 
будет иметь вид: 


ЛА + с (тА^) . (т’лХ) = М*+КМУ, (10) 
ИЛИ В неподвижной системе координат Хх: 
Х=А'/" (КРМ + М° -6 (ТА. (Т'^*.) |. (11) 


Уравнения движения обобщенной роторной системы. Обозначим точки закрепления вала в 
3 3 

подшипниках как а= » ха, иь=» хЬ (рис. 2). Кроме того, будем полагать, что а; и Ь, оста- 
1] 1=1 

ются неизменными независимо от вращения основного вала. Обозначим вектор силы, порождае- 

мый динамической связью (силовой функцией), в точках а и Ъ соответственно Ё“ и Е”, а сум- 


марный момент трения — М“. Принимая / за длину вала и используя матрицу Р, несложно по- 
лучить значения координат а и Ь: 


о ЦАМ+АА,) оо КАЖ +АЬ) 
0 1 2 3 0 1 2 3 
ЦА. Л. — АА, ДА. Л. — № А, 
ее т т 
а, =с, +1/2 Ь, =с. —1/2. 





Е’ (а,а',с) у Е“ (а,а",с) 







МЫ, с) М’ 
а! М“«а,а',с) 







ЕЬ,Ъ',с) 


| хз. 
#9 


Е ®,Ъ,с) я | аз=2 


Рис. 2. Обобщенная роторная система 


Будем считать, что на основной вал действует суммарный, внешний по отношению к ро- 


торной системе, момент М’ и сила Е’. Тогда, учитывая, что суммарный момент, действующий 
на основной вал, можно определить по формуле: 


М =М’ +М* +ах Е‘ +ьхР?, (13) 
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и используя (12), получаем: 
МХ = Е’ - Е“ )+М'+М®, (14) 
де М’, М®, Е’, Е" - вектор-столбцы с координатами соответствующих векторов в Х; 
[, — матрица длин, определяемая как 
21(^../.. РА ) 


0 ЕВ 
Е т 
= _ й 21(. +2.) 


2—2 —12 +22 
2Ц^.^. — №№ ) 21(^.^. + №1. ) 
222—242 ААА 
Исходя из вышесказанного получим систему уравнений, определяющих динамику ротор- 
ной системы: 


0 


АР. (КРЕ(Е - *)+м* чм’ |+м'-св (ТА) .(т” |}, 
(15) 
в=— (Е +Е-Е’). 
т 
Для окончательного построения модели роторной системы, как видно из (15), необходимо 


задать РЁ” ‚Е, М — М’, Е’. И если последние два параметра являются произвольными 
функциями, определяющими некоторое внешнее воздействие на систему, то первые три парамет- 
ра определяются динамической связью, входящей в состав роторной системы. Очевидно, что они 


являются функциями состояния роторной системы /, и с.. Задание этих функций определит 


окончательный вид математической модели роторной системы. В частности, Е”, Е“, М“ можно 
определить, основываясь на рассмотрении подшипниковых креплений ротора как трибосистем, по 
аналогии с тем, как это сделано в исследовании [3]. 

Приведем в качестве примера результат моделирования (рис. 3). При этом рассматрива- 
лось, что ротор закреплен в низкоскоростных подшипниках скольжения, а в качестве внешнего 
источника момента был использован двигатель постоянного тока. Качественный вид сечений ро- 
тора и статора в контактах а и Ъ приведен на рис. 4 (в контакте а наблюдается импульсный де- 
фект ротора и статора). 
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-008 008 054 002 0 002 0404 006 008 01 а), мкм 270 


а) 6) 
Рис. 3. Результаты моделирования: а — движение ротора в плоскости а,а, в контакте а; 
6 — зависимость частоты вращения ротора от его угла поворота 
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Контакт а Контакт Б 


Рис. 4. Качественный вид ротора и статора в подшипниковых узлах 


Из рис. 3 видно, что наличие импульсного дефекта оказывает существенное влияние не 
только на движение ротора в плоскости контакта а , но и на угловую частоту вращения ротора. Кроме 
того, виден периодический характер влияния дефекта в зависимости от угла поворота вала «. 
Выводы. Предложенная формула Эйлера в параметрах Родрига-Гамильтона, позволяет осущест- 
влять математическое моделирование пространственного движения твердых тел. При этом рас- 
смотрение осуществляется в неподвижном базисе, что облегчает использование данной формулы. 
Кроме того, с точки зрения численного моделирования, ее несомненным преимуществом является 
отсутствие тригонометрических функции. 

На основе формулы Эйлера в параметрах Родрига-Гамильтона построена модель ротор- 
ной системы, позволяющая анализировать ее динамические свойства. Подобная модель может 
быть использована для решения различных задач, в частности задач технической диагностики. 
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